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NOMENCLATURE 


Slope  for  vs.  7^.  Taken  here  as  2.5  (1/0.4). 
y-intercept  for  vs.  y* 

Chord  length 

5* 

Constant  for  power-law  curve-fit,  —  =  C  x  Re^ 

c 

Drag 

Drag  coefficient  of  section, - ^ — 


Portion  of  drag  coefficient  due  to  pressure  drag 
Portion  of  drag  coefficient  due  to  fiiction  drag 


Skin  fnction  coefficient. 


Local  skin  friction  coefficient 
jy-intercept  for  u  vs.  In  y 

5' 

Shape  factor,  =  — 


Hydraulic  smoothness,  = 


Surface  roughness  (in  meters) 

Arbitrary  characteristic  length  scale 

Slope  of  logarithmic-linear  relationship  u  vs.  Iny 

5* 

Power  for  power-law  curve-fit,  —  =  C  x  Re~^  " 

c 

Uc 

Reynolds  number  based  on  chord,  =  — 


Reynolds  number,  = — ,  where  /  can  be  either  c,  x,  or  6. 
V 

Mean  axial  velocity  u(y) 


Friction  velocity,  =  ,/ - 

V  P 
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Non-dimensional  mean  velocity  in  inner  variables,  =  — r 

u 

Reference  velocity 

Velocity  at  edge  of  boundary  layer,  where  y  =  8 

Distance  from  leading  edge  of  foil  to  measurement  location. 

Distance  from  wall  (in  meters) 
yu* 

Inner  variable,  = - 
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Kinematic  viscosity 
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ABSTRACT 

Velocity  profiles  were  measured  on  a  large  hydrofoil  model  to  investigate  the 
Reynolds  number  dependence  in  a  range  of  6x10^  to  6x1 0^.  The  friction  velocity 
and  skin  friction  are  computed  directly  from  the  velocity  measurements  using  the 
universal  velocity-distribution  law.  The  boundary  layer  profiles  are  also 
numerically  integrated  to  yield  displacement  and  momentum  thickness,  as  a 
function  of  Reynolds  number.  Measurements  from  the  foil  pressure-side,  which 
is  essentially  a  flat  plate,  are  shown  to  be  in  good  agreement  with  skin  friction 
measurements  at  the  same  Reynolds  numbers.  A  scaling  analysis  of  the 
measurements  reveals  approximately  a  l/l(f^  power-law  relationship  of  the 
displacement  thickness  with  Reynolds  number. 

ADMINISTRATIVE  INFORMATION 

The  Office  of  Naval  Research  funded  the  University  of  Michigan  to  perform  the 
experiments  under  contract  numbers  NOOO 14-99-1 ’0341  and  NOOO 14-99- 1-0856.  Jude  Brown,  of 
Code  508,  sponsored  the  Reynolds  number  scaling  analysis  of  the  data  presented  in  this  report. 
The  scaling  analysis  was  conducted  by  the  Naval  Surface  Warfare  Center,  Carderock  Division 
(NSWCCD),  Hydromechanics  Directorate,  Propulsion  and  Fluid  Systems  Department  (Code 
5400)  under  work  unit  number  01 -1-5080-209-32. 


6 


NSWCCD-50-TR-2001/064  August2001 


INTRODUCTION 

In  response  to  an  Office  of  Naval  Research  (ONR)  request  to  address  basic  research  in 
the  area  of  High  Reynolds  number  flows,  a  three-phase  experimental  effort  was  undertaken  at  the 
Large  Cavitation  Channel  by  the  University  of  Michigan.  The  overall  goal  of  this  experimental 
effort  was  to  increase  our  knowledge  of  turbulence  in  high  Reynolds  number  flows  [1].  The 
results  of  this  research  will  help  to  improve  prediction,  analysis,  and  control  of  full-scale 
turbulent  flows  about  ships  and  submarines. 

The  results  presented  herein  seek  to  reveal  findings  from  the  phase  I  experiment, 
specifically  relating  to  the  scaling  of  model-scale  to  full-scale  boundary  layers.  Not  until  recently, 
have  facilities  been  capable  of  testing  near  full-scale  Reynolds  number  ranges.  Therefore,  little  to 
no  experimental  data  exist  for  boundary  layer  measurements  at  Reynolds  numbers  in  the  10^ 
range.  This  report  reveals  the  Reynolds  number  relationship  to  the  boundary  layer  at  chord 
Reynolds  numbers  range  from  6x  1 0®  to  6x  1 0’. 

Blasius’  1/5*  power-law  relationship  for  displacement  and  momentum  thickness 
(deduced  from  the  l/T^-power  velocity  distribution  law)  only  gives  good  agreement  with 
experimental  data  up  to  Reynolds  numbers  of  10*.  This  report  shows  experimental  evidence  for  a 
1/10*  power-law  relationship  for  the  tested  Reynolds  number  range. 

EXPERIMENT 

Experiments  were  performed  at  the  William  B.  Morgan  Large  Cavitation  Channel  (LCC) 
in  Memphis,  TN.  As  illustrated  in  Figure  1,  the  LCC  is  a  closed-loop  re-circulating  water 
channel.  It  has  a  3-meter  square  test  section,  and  is  about  13  meters  long.  The  channel  has  a 
speed  range  from  0.5  to  18  m/s.  The  ambient  pressure  in  the  test  section  can  be  varied  from  3.5 
to  414  kPa.  Etter  and  Wilson  [2]  have  discussed  the  Large  Cavitation  Channel  in  detail. 

Moreover,  Blanton  [3]  describes  measurement  capabilities  and  provides  an  uncertainty  analysis  of 
velocity  and  pressure  measurements  taken  at  the  LCC. 

The  test  device  was  a  foil,  named  “HIFOIL”  (High  Reynolds  Number  FoU),  which  is 
typical  of  a  full-scale  section  of  a  propeller.  The  foil  is  a  modification  to  a  NACA  16  thickness 
form  with  a  =  0.8  mean-line  camber.  The  aft  3%  of  the  suction-side  is  modified  to  intentionally 
separate  the  flow,  giving  the  foil  an  “anti-singing”  characteristic.  The  pressure  side  of  the  foil  is 
completely  flat  aft  ofx/c  =  0.28.  Figure  2  shows  the  geometry  of  the  foil,  and  Figure  3  shows  the 
trailing  edge  in  detail.  HIFOIL  has  a  chord  length  of  2.134  meters  and  a  maximum  thickness  of 
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0.171  meters.  The  surface  roughness  is  less  than  0.25 fm,  which  corresponds  to  a  maximum 
hydraulic  smoothness  k*  =  0.2  at  18.3  m/s.  At  the  highest  speed,  18.3  m/s,  the  foil  generated  an 
estimated  lift  force  of  670  kN  (150,000  Ibf)  [4].  The  foil  was  placed  at  a  zero  angle  of  attack. 

Laser  Doppler  Velocimetry  (LDV)  measurements  were  taken  on  the  pressure  side  at  four 
xtc  locations:  0.930, 0.957,  0.979,  and  1.0.  On  the  suction  side,  stations  were  analyzed  at  three 
xlc  locations:  0.930, 0.957,  and  0.979.  In  addition,  a  station  in  the  wake  atx/c  =  1.028  was  also 
investigated.  The  experiment  was  performed  at  four  different  speeds:  3,  6, 12,  and  18  m/s.  These 
correspond  to  nominal'  Reynolds  numbers  of9.0xl0^  1.9xl0’,  3.7xl0’,  and  5.6xl0’ respectively. 
The  foil  was  allowed  to  naturally  transition  fi-om  laminar  to  turbulent  flow.  This  fact  has 
introduced  some  small  scatter  in  the  data  shown,  because  the  transition  location  probably  changes 
with  speed.  This  scatter  is  most  evident  in  the  3  and  6  m/s  data.  Future  experiments  plan  to 
incorporate  trip-wires  to  force  transition  on  the  foil  at  a  fixed  location.  Bourgoyne  [4]  describes 
this  experiment  in  detail  and  presents  turbulence  measurements  for  the  current  data. 


SCALING  ANALYSIS 

The  pressure-side  of  the  foil  is  assumed  to  have  the  same  flow  physics  as  flat-plate  flow. 
This  hypothesis  was  verified  three  ways:  (1)  verifying  the  pressure  gradient  is  zero  on  the 
pressure-side  of  the  foil,  (2)  comparing  local  skin  fiiction  coefficients  with  flat-plate  theory,  and 
(3)  verifying  that  profiles  taken  at  different  x/c  locations  exhibit  shape  factor  similarity. 
Bourgoyne  [4]  measured  the  pressure  distribution  on  both  the  pressure  and  suction  side  of 
HIFOIL  during  phase  II  of  the  experiment.  The  pressure-side  results  do  indeed  show  a  negligible 
pressure  gradient  in  the  trailing  edge  region,  where  the  velocity  profiles  are  measured.  Local  skin 
friction  and  shape  factor  correlations  will  be  given  in  the  Skin  Friction  and  Characteristic 
Numbers  sections,  respectively. 

To  investigate  the  Reynolds  number  dependence  of  the  boundary  layer,  the  velocity 
profiles  are  numerically  integrated  to  compute  displacement  and  momentum  thickness  and  the 
shape  factor.  Also,  the  profiles  are  plotted  in  inner-variable  coordinates,  vs.  which  clearly 
show  the  effect  of  Reynolds  number  on  the  profiles.  The  skin  fnction  coefficients  presented  are 
not  measured,  but  instead  are  inferred  from  the  velocity  measurements.  The  Logarithmic-linear 
Regression  section  explains  how  this  analysis  was  done. 


'  The  tenqjerature  of  throughout  the  experiment  varied  from  22  to  40  degrees  Celsius  (71.5  to  104.2 
degrees  Fahrenheit),  causing  the  Reynolds  number  to  vary  for  a  fixed  speed.  Although  for  the  data  given 
here,  the  Reynolds  number  remained  relatively  constant  for  each  condition,  and  should  be  taken  as  those 
given  above. 
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VELOCITY  PROFILES 


As  shown  in  Figure  4,  the  HIFOIL  boundary  layer  profiles  show  negligible  variation  with 
Reynolds  number,  when  looking  at  mean  velocity  profiles.  Normalizing  the  profiles  with  the 
edge  velocity,  Ue,  instead  of  C/*,  reveals  a  small  variation  in  the  boimdary  layer  region  and 
highlights  that  the  changes  are  occurring  in  the  overlap  layer.  For  example,  compare  Figure  4d 
and  Figure  5.  The  changes  become  obvious  when  the  data  are  normalized  by  the  fiiction  velocity, 
u*.  Figure  6  shows  the  experimental  measurements  from  the  pressure-side  plotted  in  inner 
variable  coordinates.  The  next  section  shows  the  method  used  to  shift  and  scale  the  boundary 
layer  data,  such  that  it  can  be  plotted  in  u*.  y*  coordinates.  In  all  the  velocity  profile  figures, 
Spalding’s  [5]  law  of  the  wall  model  is  shown.  Spalding’s  single  formula  for  the  law  of  the  wall 
[5]  is  given  as: 


y  =u  +e 


■kB 


-l-KU* 


2  6 


(1) 


where  k  is  a  dimensionless  constant,  taken  here  as  0.40.  u*  andy^  are  inner-variable  coordinates. 
All  profiles  collapse  to  Spalding’s  law  of  the  wall  model,  though  each  profile  was  shifted 
independently. 

Figure  7  shows  axial  mean  velocity  profiles  in  outer  variables  for  the  suction-side  of  the 
foil,  taken  at  x/c  =  0.930,  0.958,  and  0.979,  respectively.  Figure  8  shows  the  same  data  plotted  in 
inner-variable  coordinates.  This  data  was  more  difficult  to  shift  correctly,  because  only  a  few 
points  were  in  the  logarithmic  region.  On  the  suction-side,  separation  appears  to  occur  just 
downstream  of  the  x/c  =  0.979  profile.  The  adverse  pressure  gradient  tends  to  shift  the  data  above 
the  logarithmic  law  [8].  In  the  x/c  =  0.979  profile  (Figure  8e),  the  data  has  not  been  shifted  due  to 
the  difficulty  in  determining  the  logarithmie  region.  Note  that  the  x/c  =  0.979  profile  was  just 
downtream  of  the  “knuekle”,  as  shown  in  Figure  3. 

LOGARITHMIC-LINEAR  REGRESSION 

This  section  explains  how  the  data  is  shifted  to  fit  the  universal  velocity-distribution  law, 
and  how  u*  and  C/are  computed.  The  logarithmic  velocity-distribution  law  states  a  linear- 
logarithmic  relationship  between  the  velocity  and  the  distance  fi’om  the  wall  as  such: 

u*=A\ny*+B  (2) 
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where  typically  A  =  2.5  and  B  =  5.5.  Here, 


M^=4and/=^whereM 

U  V  V  p 


M*  is  the  friction  velocity,  and  is  a  measure  of  the  turbulent  edd)dng  in  the  outer  boundary  layer 
[6].  Also,  To  is  the  shear  stress  at  the  wall  and  p  is  the  density  of  the  fluid.  Thus,  using  the 
previous  equations,  a  new  dimensional  form  of  the  universal  velocity-distribution  law,  is: 


M  =  -\rBu 

\l  J  I  V  j 


Let  y  =  z  —  Zg,  where  Zg  represents  the  location  of  the  wall  surface.  Since  the  precise  location  of 
the  wall  surface  is  not  known,  it  must  be  determined  iteratively.  Furthermore,  define 

f  lu^ 

m  =  Au  and  d  =  Au*  In  —  +Bu’ as  shown  above,  which  leaves 

I  V  j 


M  =  w  In  - — —  +  d 

I  /  J 


Here  u  and  z  are  given  (in  m/s  and  m  respectively)  from  the  experiment,  /  is  an  arbitrary  scale 
length  (taken  here  as  one),  and  m  and  d  are  determined  from  a  linear  regression  analysis  of  the 
linear-log  region  of  the  boundary  layer.  The  Zo  is  determined  iteratively  such  that  B=5.5,  now 
given  by 

B  =  ^-A\n—  (6) 

U  V 

The  friction  velocity  can  now  conveniently  be  computed  from: 


*  m 

u  =  — 
A 


Also,  mdy^  are  computed  from  the  friction  velocity,  w*,  in  addition  to  the  local  skin  friction 

coefficient,  which  is  computed  as: 

c.=Tj(pU^/2)  =  (u/Uy  (8) 
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Table  1  shows  the  distance  in  mm  that  each  velocity  profile  was  shifted.  For  most  cases,  the  shift 
was  smaller  than  the  probe  diameter  (0.2mm)  of  the  LDV. 


Table  1.  Distance  (in  mm)  that  data  is  shifted  such  that  B  =  5.5 


x/c 

6 

12 

18 

Pressure 

0.930 

0.37 

0.28 

0.18 

0.47 

0.958 

0.38 

0.08 

0.17 

0.30 

0.979 

0.40 

0.03 

0.05 

0.17 

1.000 

-0.09 

-0.03 

-0.08 

-0.07 

Suction 

0.930 

-0.15 

0.11 

-0.02 

-0.31 

0.958 

-0.05 

-0.23 

-0.29 

0.00 

SKIN  FRICTION 

Various  skin-fnction  formulas  are  given  in  Schlichting  [6],  each  tending  to  work  well  for 
a  specified  Reynolds  number  range.  Three  are  applicable  to  Reynolds  number  range  of  the 


HIFOIL  experiment: 

Blasius  Theory 

=  0.0592(i?eJ-'/' 

5x10*  <  Re,  <  10’ 

(9) 

Prandtl-Schlichting 

=21og(i?e,-0.65)"'^ 

10’ <  Re,  <10® 

(10) 

Schultz-Grunow 

=  0.0592(i?eJ"'/^ 

10*  <  Re,  <  10* 

(11) 

Figure  9  shows  the  local  skin  friction  coefficient  from  the  foil  pressure-side  compared  with  the 
previous  formulas.  The  HIFOIL  experimental  data  is  shown  to  lie  below  the  theoretical  models. 
This  is  most  likely  due  to  a  difference  in  transition  locations  between  experiment  and  theory.  In 
fact,  the  Prandtl-Schlichting  formula  (for  total  skin  fnction)  even  includes  a  correction  for 
transition  location.  Kempf  [6]  achieved  experimental  total-skin-fiiction  coefficients  up  to  Re  = 
5x10*.  His  data  tended  to  lie  between  midway  between  Equations  (10)  and  (1 1)  for  the  Reynolds 
number  range  in  this  experiment. 

CHARACTERISTIC  NUMBERS 

To  show  Reynolds  number  effects,  the  displacement  thickness,  momentum  thickness  and 
shape  factors  are  computed. 
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Displacement  Thickness 


The  displacement  thickness  is  a  characteristic  number  representing  the  distance ' 
streamlines  are  displaced  due  to  the  presence  of  the  boundary  layer,  and  is  defined  as: 


It  was  integrated  numerically  using  the  Trapezoidal  rule.  Experimental  measurements  taken  with 
the  LDV  yielded  average  minimum  values^  of  about  80.  To  accurately  compute  the 
displacement  thickness,  Spalding’s  wall  model.  Equation  (1),  was  used  for  0  <y^<  80.  Thus, 


Spalding  Experiment 


where  y„,v,  corresponds  to  the  distance  of  the  measurement  closest  to  the  wall. 


Momentum  Thickness 


The  momentum  thickness,  like  the  displacement  thickness,  is  a  characteristic  number  that 
gives  a  measure  of  the  loss  of  momentum  in  the  boundary  layer,  owing  to  the  frictional  losses 
caused  by  the  fluid  interaction  with  the  wall.  A  similar  computation  is  performed  to  compute  the 
momentum  thickness,  defined  as: 


0=  J  1- 


u  u 

u^iTT. 


u 

\  “ 

Ue, 

Similarly,  this  can  be  written  as 


Spalding  Experiment 


(15) 


^  The  distance  to  the  wall  was  limited  to  the  accuracy  of  the  laser.  The  LDV  has  a  probe  volume  of  02mm 
which  corresponds  to  wall  coordinates  between  60  and  160. 
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Figure  10a  shows  the  results  of  this  research  for  the  pressure-side  data.  Notice  the  data 
scatter  for  the  four  different  Reynolds  numbers.  The  3  m/s  (Re  =  9x10*)  shows  little  data  scatter, 
and  the  resulting  displacement  thickness  is  less  than  the  other  cases.  The  transition  location 
possibly  is  further  aft  in  the  lowest  Reynolds  number  case,  causing  “thinner”  results  for  the  lower 
Reynolds  number  case.  The  6  m/s  data  (Re  =  19x10*)  exhibits  more  data  scatter,  which  can  be 
attributed  to  transition  location  uncertainty.  On  the  other  hand,  the  high  Reynolds  number  case, 

18  m/s  (Re  =  5.6xl0’)  shows  to  be  fully  turbulent. 

Figure  10b  indicates  the  suction-side  displacement  and  momentum  thickness  as  a 
function  of  Reynolds  number.  The  data  does  not  collapse  as  in  Figure  10a,  due  to  the  changing 
pressure  gradient  in  the  trailing  edge  region. 

The  Reynolds  number  is  computed  based  on  the  momentum  thickness.  This  is  defined 
as: 

u  e 

(16) 

V 

Figure  1 1  shows  a  log-log  relationship  between  Ree  and  Re^. 

A  power-law  curve  fit  of  the  displacement  thickness  versus  Reynolds  number  for  the 
pressure-side  data  is  shown  in  Table  2.  The  average  computed  in  Table  2  does  not  include  the 
data  in  the  wake  {x/c  =  1 .028).  Here,  C  and  n  are  the  constants  in  the  following  expression: 

—  =  CxRe~^"'  (17) 


c 

Table  2.  Power-law  constants  for  computing  <5*  on  foil  pressure-side. 


x/c 

c 

n 

Blasius 

0.046 

5 

0.930 

0.0088 

10.6 

0.956 

0.0106 

9.5 

0.978 

0.0119 

8.9 

1.000 

0.0109 

9.2 

1.028 

0.0272 

9.8 

Average  (excluding  x/c  =  1.028) 

0.0106 

9.6 
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Shape  Factor 


The  shape  factor  is  a  dimensionless  characteristic  number  quantifying  the  velocity  profile 
shape,  and  is  given  here  as: 

H=—  (18) 

0 

The  shape  factor  is  heavily  dependent  on  the  external  pressure  gradient  [6].  Hinze  [7]  reports  a 
dependence  of  shape  factor  on  Reynolds  number,  but  the  variation  is  small  for  high  Reynolds 
numbers.  Based  on  Blasius’  theory  [6],  the  shape  factor  for  the  flat  plate  is  independent  of 
Reynolds  number  and  has  a  value  of  1.278.  Figure  12  shows  the  shape  factor’s  dependence  on 
Reynolds  number  for  both  the  pressure  and  suction-side  experimental  data.  As  is  expected  the 
pressure-side  data  shows  little  dependence  on  the  Reynolds  number  and  collapses  on  the 
theoretical  prediction.  This  again  is  verified  in  Figure  13  where  four  different  profiles  collapse, 
each  profile  being  scaled  by  its  edge  velocity.  This  shows  that  the  pressure-side  of  HIFOIL  does 
produce  very  similar  flow  physics  to  that  of  a  flat  plate.  On  the  other  hand,  the  suction-side  data 
tends  to  change  drastically  depending  on  the  x/c  location  where  the  data  were  measured.  The  x/c 
=  0.978  shape  factor  is  much  higher  than  the  two  cuts  at  x/c  =  0.9S0  and  x/c  =  0.958,  owing  to 
the  fact  that  this  profile  is  just  downstream  of  the  knuckle,  as  shown  in  Figure  3.  This  indicates 
the  strong  dependence  on  the  external  pressure  gradient  that  was  alluded  to  earlier.  Previous 
research  suggests  separation  occurs  at  a  value  of  7/ of  about  4.0  [8]. 


SECTION  DRAG 

The  drag  on  the  hydrofoil  is  directly  related  to  the  momentum  deficit  in  its  wake.  This 
can  be  expressed  as: 

Co  =2^  (19) 

c 

where  is  the  momentum  thickness  at  a  distance  downstream  of  the  body,  computed  from 

Equation  (6).  Figure  14  shows  the  LDV  measurements  taken  at  x/c  =  1.028.  These  data  have 
been  normalized  by  the  edge  velocity.  Because  of  the  nature  of  the  separated  flow  from  the  blunt 
trailing  edge,  the  data  in  the  wake  exhibit  higher  turbulence  for  the  suction  side  of  the  foil  (y  >  0), 
and  thus  shows  more  scatter  in  the  data. 
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In  order  to  investigate  Reynolds  number  scaling  from  the  measured  drag,  the  base  drag 
must  be  subtracted  from  the  total  drag.  For  this  case,  Qwas  estimated  using: 


ic,(Re) 


0.074 


(20) 


where  the  14  accounts  for  just  one  side  of  the  foil.  Hence,  to  obtain  the  total  skin  friction  for  both 
sides  of  the  foil.  Equation  (20)  is  doubled.  Unless  the  separation  point  moves,  pressure  drag 
should  be  a  constant  value  (independent  of  Reynolds  number).  Hence,  an  average  value  of 

was  used  to  estimate  the  portion  of  drag  due  to  pressure  loss,  as  such: 

=  Co  -|(C^(1.9x10’)  +  C^(3.7a:10’)  +  C^(5.6a:10’))  (21) 

Cop°  here  is  computed  to  be  0.0054  and  is  subtracted  from  the  drag  coefficient  computed  in 
Equation  (8)  to  compute  Qy,  the  portion  of  drag  due  to  friction: 

Cof=Co-Ct7  (22) 

Figure  15  shows  the  results  of  this  study.  Most  of  the  pressure  drag  is  due  to  the  flow 
separation  on  the  suction-side  of  the  foil.  It  is  interesting  to  note,  that  when  the  pressure  drag  is 
subtracted  from  the  total  drag,  the  resulting  power,  n,  is  9.6,  which  is  almost  identical  to  the 
power  previously  shown  in  the  boundary  layer  studies.  Shen  developed  a  theory  for  more 
accurately  predicting  «,  for  a  given  Reynolds  number  ratio.  For  the  Re)molds  number  range 
presented  here,  Shen’s  theory  predicts  n  =  9.5  [9].  This  is  almost  exactly  what  the  experiment 
shows  in  Table  2. 


CONCLUSIONS 

Boundary  layer  data,  measured  in  a  Reynolds  number  range  from  6x10*  to  6xl0’,  have  shown  to 
exhibit  a  1/10*  power-law  relationship  with  Reynolds  number.  Skin  friction,  inferred  from 
measurements  taken  on  the  pressure-side  of  the  foil,  is  in  good  agreement  with  existing  data  taken 
at  the  same  Reynolds  number.  Displacement  and  momentum  thickness  are  computed  by 
numerically  integrating  the  velocity  profiles.  The  shape  factor  is  also  computed  and  shows  shape 
similarity  for  the  four  locations  on  the  pressure-side  of  the  foil.  The  momentum  thickness  in  file 
wake  is  used  to  compute  the  drag  coefficient  on  the  foil.  Frictional  drag  is  deduced  from  the  total 


NSWCCD-50-TR-2001/064  August2001 


15 


drag  and  is  shown  to  exhibit  similar  trends  as  the  measurements  taken  on  the  flat/pressure  side  of 
the  foil.  Future  experiments  should  be  performed  at  higher  Re)aiolds  number  ranges  on  flat-plate 
flow  to  show  a  higher-order  correlation  of  the  boundary  layer’s  dependence  on  Reynolds  number. 
Moreover,  more  analysis  should  be  performed  on  the  suction-side  data  to  investigate  high 
Reynolds  number  effects  of  separated  flows. 
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Figure  4.  Pressure-side  mean  axial  velocity  profiles. 


a.  x/c  =  0.930  pressure  side. 


b.  x/c  =  0.958  pressure  side. 
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Figure  6.  Pressure-side  inner-variable  velocity  profiles. 


y+ 


a.  x/c  =  0.930  pressure  side. 


b.  x/c  =  0.958  pressure  side. 
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Figure  7.  Suction-side  mean  velocity  profiles. 
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a.  x/c  =  0.930  suction  side. 
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b.  x/c  =  0.958  suction  side. 
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Figure  9.  Local  skin-friction  coefficient  computed  from  pressure-side  velocity  profiles. 
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Figure  10.  Displacement  and  momentum  thickness  versus  Reynolds  number. 


a.  Pressure-side  displacement  and  momentum  thickness  versus  Reynolds  number. 


b.  Suction-side  displacement  and  momentum  thickness  versus  Reynolds  number. 
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Figure  12.  Shape  Factor  as  a  function  of  Reynolds  number 
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Figure  13.  Pressure-side  mean  axial  velocity  profiles  showing  shape  similarity. 
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Figure  14.  LDV  wake  measurements  taken  at  x/c  =  1.028. 
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Figure  15.  Reynolds  number  dependence  of  drag  coefficient. 
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